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INTRODUCTION 


We  consider  the  case  of  a  monochromatic  cylindrical  beam  with  Gaussian  intensity 
profile,  which  is  incident  endwise  upon  a  semi-infinite  planar  waveguide  structure.  This 
study  can  have  practical  implications  for  electro-optic  devices  where  coupling  of  laser 
energy  into  guided  wave  devices  is  of  interest. 

Shown  schematically  in  Figure  1,  the  Gaussian  beam  is  incident  at  angle  0  onto  the 
z  =  0  plane  where  this  plane  defines  the  boundary  between  the  superstate  and  substrate. 
The  problem  analyzed  here  is  two-dimensional,  where  the  incident  beam  and  waveguide 
structure  are  invariant  in  the  y  dimension. 


FIGURE  1 .  Schematic  of  a  Gaussian  Beam  Incident  Endwise  Upon  a 
Semi-Infinite  Planar  Waveguide.  The  Gaussian  beam  and  waveguide 
channel  are  invariant  in  the  y-direction.  The  shaded  substrate  regions 
with  permittivity  £,  border  the  planar  waveguide  channel  which  has 

permittivity  £2-  The  substrate  regions  are  semi-infinite  in  the  ±x 
directions  and  the  waveguide  has  width  d.  Also,  the  substrate  and 
waveguide  are  both  semi-infinite  in  the  -z  direction. 
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METHOD  OF  COMPUTATION 


The  plane  z  =  0  separates  two  semi-infinite  media.  The  superstate  region  z  >  0  is 
homogeneous  and  the  substrate  region,  which  is  invariant  in  the  -z  direction,  contains  a 
planar  waveguide  structure  of  width  d.  We  assume  a  monochromatic  Gaussian  intensity 
profile  beam,  characterized  by  width  a  and  wavelength  X,  is  normally  incident  on  the 
waveguide  structure.  We  approximate  the  amplitude  of  the  Gaussian  beam  at  z  =  0  by 


P^"{x,z  =  0)  =  exTp[-x^  /o-^] 


(1) 


where  Py’'^  represents  the  y  component  of  either  the  electric  or  magnetic  field,  whichever  is 

appropriate.  When  the  (electric)  magnetic  field  is  parallel  to  the  waveguide  boundaries,  P^'^ 

is  the  (electric)  magnetic  field  amplitude.  We  wish  to  calculate  the  transmitted  and  reflected 
fields  where  we  expect  at  least  part  of  the  transmitted  energy  to  be  in  the  form  of  a  guided 
wave. 

The  approach  used  in  this  work  has  similarities  to  previous  work  (References  1  through 
5).  In  References  1  through  3,  solutions  were  obtained  as  modal  expansions  similar  to  that 
done  here.  In  Reference  1,  the  solution  was  obtained  entirely  in  Fourier  space  and 
digitization  of  real  space  is  not  done.  In  References  2  through  3,  real  space  is  discretized  as 
is  done  in  the  present  work.  In  Reference  4  and  also  in  real  space,  a  numerical  inte^ation 
of  Maxwell’s  equations  is  done  without  the  need  for  modal  expansion  solutions.  Finally, 
in  Reference  5,  the  present  method  was  applied  to  investigate  transverse  localization  of 
electromagnetic  propagation  in  alternating  high-low  dielectric  media  of  random  thickness. 
Like  the  work  here,  the  incident  beam  was  Gaussian  and  incident  endwise  on  the  laminate 
structure. 

Along  the  X  direction,  a  discrete  set  of  points,  separated  by  Ax,  is  defined  that  are 
symmetrically  located  about  x  =  0.  This  is  shown  in  Figure  2.  The  total  width  of  these 
points  is  L  =  n^Ax  and  this  defines  the  lateral  limits  (x  =  ±  L/2)  of  the  calculation  region 
along  the  x  direction.  Because  of  this,  we  require  that  L  »  a  so  that  the  incident  field  is 

essentially  zero  at  the  limits  x  =  ±  L/2.  No  attempt  is  made  to  incorporate  transparent 
boundary  conditions  at  the  lateral  limits  of  the  calculation  region.  The  field  at  the  lateral 
limit  and  beyond  is  always  assumed  to  be  zero.  If  nonzero  field  intensity  reaches  the  lateral 
limit  boundary,  for  whatever  reason,  then  the  calculation  is  terminated.  A  simple  example 
is  when  the  substrate  region  is  a  homogeneous  dielectric  material  (no  waveguide  channel). 
In  this  case,  the  width  of  the  transmitted  field,  relative  to  the  initial  Gaussian  width,  would 
spread  out  with  increasing  depth  by  means  of  diffraction  and  eventually  reach  the  lateral 
limit  boundary.  In  this  case,  the  field  would  probably  “reflect”  from  the  boundary,  and  this 
would  be  nonphysical  since  the  boundary  is  fictitious. 
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PTGURE  2.  Schematic  Illustrating  x  Coordinate  Discretization  (Solid  Dots 
Separated  by  Ax)  With  x-Limits  on  Extent  of  Calculation  Region  (Shown  by 
Length  L  Between  the  Vertical  Dotted  Lines).  The  incident  field  is  negligible  at  the 
extreme  limits.  In  this  example  the  solid  dots  represent  n^=21  digitized  x 
coordinates. 


Solutions  to  Maxwell’s  equations  in  the  homogeneous  z  >  0  region  are  straightforward 
and  not  discussed  here.  In  the  inhomogeneous  z  <  0  region,  the  solutions  are  more 
complicated  and  with  the  x  coordinate  digitized,  we  write  the  x  dependence  of  the  two 
Maxwell’s  equations  VxE(r)  =  i(C[)/c)B(r)  and  VxB(r)  =  -i(£()/c)D(r),  where 

D(r)  =  £(r)E(r),  in  centered  finite-difference  form.  In  the  z  <  0  region,'this  yields: 


dz  c  ^ 


ic 


CO{Axr 


By{x  +  Ax,z)-  By{x,z)  ^  Ry(x- Ax,z)--By(x,z) 


e(x-i- Ax/2) 


e(x- Ax/2) 


(2a) 


dEy{x,z) 

dz 


-—BJx,z) 


c 


(2b) 
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^L^  =  -1^S(x)E/x,z) 
dz  c  ^ 

^  (oilxf^Ey^^'Z)  -  Ey(x  +  Ax,z)  -  Ey(x  -  Ax,z)}  (2c) 


dHy(x,z) 

dz 


—s{x)E,(x,z) 

c 


(2d) 


where  the  z  coordinate  remains  a  continuous  variable.  The  z  components  of  the  fields  have 
been  eliminated  since  the  z  <  0  region  is  invariant  along  the  z  direction.  The  permittivity 
£(;c)  describes  the  substrate  region  where  s(x)  =  ei  when  \x\>d/2  and  £(x)  =  e2  when 
\x\<d/2.  Equation  2  is  versatile  since  the  only  part  that  depends  on  the  waveguide 
structure  is  the  permittivity  £(x),  and  this  quantity  can  be  configured  to  represent  a  variety 
of  situations.  Equation  2  is  valid  at  a  discrete  point  x  in  the  set  of  points  and  since  we 
must  collectively  consider  all  of  the  jc  points,  this  results  in  a  coupled  equation  system 
totaling  4«^  equations.  The  equations  associated  with  all  points  which  cover  distance 
L  =  n^Ax  may  be  concisely  written  in  matrix  form  as 


^'^(X’^)-  =  M(X)A(X,z) 
dz 


(3a) 


where 


fE,(X,z)A 


A(X,z)  = 


EyiX,z) 

B,(X,z) 

^By(X,z)y 


(3b) 


Since  M(X)  is  independent  of  z,  the  solution  for  the  fields  in  the  z  <  0  region  is 
straightforward  by  diagonalization  of  M(X): 


A(X,z)  =  S(X)exp(Az)C  = 


0  Yc.l 


vC  y 


(4) 
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The  square  matrix  S(X)  has  columns  which  are  the  eigenvectors  of  M(X),  exp(Az)  is  a 

diagonal  matrix  with  A  representing  the  set  of  eigenvalues  associated  with  M(X),  and  C  is 
a  column  vector  of  constants.  In  the  right  hand  side  of  this  equation,  the  S  and  exp 
matrices  have  been  sectored.  The  exp  matrix  has  exponential  terms  along  the  diagonal 

where  the  set  of  eigenvalues  A  consists  of  pairs  such  that  for  every  eigenvalue  +X,  there  is 
another  eigenvalue  which  is  equal  but  opposite  in  sign,  i.e.  -X.  The  eigenvalue  pairs  +?i 

and  -X  are  associated  with  solutions  that  propagate  or  decay  exponentially  in  the  -z  and  +z 
directions,  respectively.  With  this  in  mind,  it  is  clear  in  Equation  4  that  the  exp  matrix  has 
been  sectored  by  separating  the  solutions,  which  propagate  into  the  ±z  directions.  For  the 
problem  considered  in  this  work,  solutions  in  the  z  <  0  region  having  coefficient  C_  are 
unphysical  since  they  represent  propagation  in  the  +z  direction  or  exponentially  increasing 
behavior  as  z  ^  -<» .  Because  of  this  we  set  C_  =  0  and  the  z  ^  0  solution  may  now  be 
written 


^E'(X,z)^ 

,B'(X,z); 


'S„(X)e^C,^ 

ts,,(X)e^cJ 


(5) 


where  we  have  introduced  the  superscript  t  to  denote  transmitted  field  and  the  tilde  notation 
to  denote  an  abbreviated  column  vector  as 


and  B  = 

[^y 

K^yJ 

(6) 


From  Equation  5  we  easily  relate  the  magnetic  and  electric  fields  in  the  z  ^  0  region  as 


B'(X,z)  =  S3,(X)[S„(X)]-'E'(X,z) 


(7) 


We  now  apply  the  boundary  conditions  at  the  z  =  0  boundary  with  the  requirement  that 
the  tangential  components  of  the  electric  and  magnetic  fields  be  continuous  across  this 

boundary.  These  conditions  can  be  written  as  E'"‘'(X,0)  +  E''(X,0)  =  E'(X,0)  and 

B'"‘'(X,0)  +  B'’(X,0)  =  B'(X,0),  where  the  superscripts  inc,  r,  and  t  denote  the  incident, 
reflected,  and  transmitted  field,  respectively.  Using  Equation  7,  these  continuity  conditions 
may  be  written  as 


^E'""(X,0)^ 

-U 

^E''(X,0)^ 

^E'(X,0)^ 

^  b:'(X,o)  ' 

,B'"^(X,0); 

,B'(X,0); 

,S„(X)[S„(X)]-'e'(X,0)^ 
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The  incident  fields  represent  the  Gaussian  incident  beam.  To  solve  for  the  transmitted  E' 

and  reflected  E'^  electric  fields,  we  need  a  relationship  between  E'’  and  B'"  and  since  these 
latter  fields  are  in  the  homogeneous  superstate  region,  we  can  write  these  fields  as  Fourier 
expansions 


E''(X,z)  =  F(X,K)E''(K,2)  and  r(X,z)  =  F(X,K)B''(K,z)  (9) 


These  Fourier  expansions  are  written  in  matrix  form  where  the  square  matrix  F(X,K)  is  a 
Fourier  transform  operator.  The  K  denotes  the  set  of  wave  vectors  in  Fourier  space 
analogous  to  the  set  of  x  coordinates  X  in  real  space.  In  Fourier  space  and  with  Maxwell’s 

equations,  we  can  further  relate  E'^(K,0)  and  B'^(K,0)  by 


B''(K,0)  =  Z(K)r(K,0) 


(10) 


where  Z(K)  is  a  square  matrix  defined  elsewhere  (References  1  and  3).  Using  Equations  9 
and  10  in  Equation  8  yields 

I  Y  E'(X,0) 

AiSn  -FZF-‘Ye'-(X,0)^ 

where  I  is  the  identity  matrix.  Equation  11  can  be  solved  for  E'(X,0)  and  then,  from 
Equation  5  we  find  that  =  S^/E'(X,0)  that  further  yields  the  transmitted  electric  field  for 
depths  z  <  0  as 


'  E'""(X,0) 

fr""rx.0'> 


(11) 


E'(X,z)  =  S„eX/E^(X,0) 


(12) 


This  is  our  primary  theoretical  result  from  which  we  obtain  numerical  results  below. 


NUMERICAL  RESULTS 


The  numerical  results  were  obtained  with  the  following  physical  parameters:  co/c  =1.5 

inverse  micrometers  (wavelength  =  4.19  micrometers),  c  =  3.0  micrometers,  d  =  3.0 
micrometers  and  L  =  30  micrometers.  The  number  of  digitized  points  was  =  1001, 
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which  yields  Ax  =  0.03  micrometers.  The  permittivity  values  are  £j  =  (-16,0)  and 
£2  =(2.25,0). 

In  Figure  3,  we  illustrate  the  existence  of  guided  wave  modes  by  looking  at  the  inverse 
square  of  the  dispersion  relation  for  both  p-  and  s-polarization.  The  dispersion  relations  are 
obtained  for  the  z  <  0  region  by  matching  boundary  conditions  and  finding  the 
electromagnetic  field  solutions  within  this  region.  The  peaks  indicate  the  wave  numbers 

which  satisfy  the  dispersion  relation  for  the  fixed  co  and  therefore  represent  wave  numbers 

associated  with  guided  mode  solutions.  For  frequency  co/c  =  1.5  inverse  micrometers  and 
d  =  3.0  micrometers,  we  see  that  there  are  three  p-polarization  mode  peaks  and  two 
s-polarization  mode  peaks.  These  peaks  indicate  guided-mode  solutions,  which  propagate 
in  the  z  direction  along  the  guided  wave  stracture  described  in  this  work.  Figures  4 
through  6  describe  coupling  of  the  incident  beam  into  these  guided  modes. 


Normalized  Wave  Number,  A:^/(co/c) 

FIGURE  3.  Guided-Wave  Resonance  Versus  Wave  Number 

at  Fixed  Frequency  ®.  This  graph  indicates  three  p-pol 
(magnetic  vector  parallel  to  guide  boundaries)  mode  solutions 
and  two  s-pol  (electric  vector  parallel  to  guide  boundaries)  mode 
solutions. 


For  p-polarization,  Figure  4  shows  the  intensity  of  the  incident  field  (at  z  =  0)  and 

transmitted  field  (Equation  12  at  z  =  0  and  z  =  -500X,).  It  is  seen  that  there  is  significant 
coupling  of  the  incident  beam  into  the  planar  waveguide  where  the  reflectance  and 
transmittance  is  0.335  and  0.665,  respectively. 
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In  Figure  4,  it  is  seen  that  there  are  some  interesting  observations  about  the  three  p-pol 
modes  shown  in  Figure  3.  First,  the  mode  indicated  at  /(co/ c)  =  0.9  (k^=l .35  inverse 
microns)  is  associated  with  a  guided  mode,  which  has  peak  intensity  around  x  =  0  and 
minimum  intensity  at  the  guide  boundaries.  This  intensity  peak  is  seen  as  the  slight 
negative  curvature  (about  jr  =  0  )  of  the  two  guided-mode  curves.  Also,  this  mode  is 
propagating  (non-evanescent)  in  both  the  x  and  z  directions.  The  other  two  p-pol  peaks  in 
Figure  3  are  centered  about  k^  /(to/c)  =  1.62 .  This  indicates  that  they  are  associated  with 
coupled  surface  waves  that  are  evanescent  relative  to  the  metal-dielectric  waveguide 
boundaries.  Thus,  these  modes  are  evanescent  in  the  x  direction  and  propagating  in  the  z 
direction.  To  see  this,  we  note  that  a  surface  plasmon  that  propagates  ^ong  a  metal- 

dielectric  interface  has  dispersion  relation  k^  /(co  /c)  =  -^£,£2  +  ^2)  for 

permittivity  parameters  used  here,  this  yields  k^p/(ty/c)  =  1.618.  The  two  p-pol  peaks 

seen  in  Figure  3  are  split  about  this  value  with  the  splitting  due  to  coupling  between  the 
surface  waves  at  each  interface.  Finally,  we  note  that  since  these  modes  are  evanescent, 
their  intensity  is  a  maximum  at  the  guide  interfaces. 


x-Coordinate 


FIGURE  4.  Electric  Field  Intensity  Versus  x  Coordinate  For 
P-Polarization  and  oo/c  =1.5  Inverse  Microns.  For  reference, 
the  z  =  0  Gaussian  beam  envelope  is  indicated  (dash  line).  The 
remaining  two  intensity  curves  are  for  the  transmitted  field  (Eq. 

(12))  at  z  =  0  (thick  line)  and  z  =  -500?i.  (thin  line).  The  planar 
waveguide  boundaries  are  at  x  =  ±  1.5  micrometers. 


For  p-polarization,  there  is  another  surface  mode,  not  indicated  in  Figure  3,  with  which 
there  is  coupling:  the  surface  plasmon  at  the  z  =  0  metal-vacuum  interface.  This  surface 
plasma  mode  is  confined  to  the  metal-vacuum  interface  (evanescent  in  the  z  direction)  and 
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propagates  in  the  x  direction.  Coupling  can  occur  between  the  incident  and  metal-vacuum 
surface  plasmon  field  because  of  the  discontinuity  of  the  dielectric  waveguide  channel. 
This  surface  plasmon,  with  wave  number  has  dispersion  relation 

k^p  l{(Olc)  =  ^£2 /(e2  + 1)  =  1.0324  for  £2  =(-16,0).  This  surface  plasmon  coupling  is 
evident  in  Fig.5  which  shows  the  Fourier  spectrum  of  the  electric  field  at  z  =  0  and  -500X. 
Figure  5  is  a  plot  of  the  Fourier  transform  of  E'(X,z),  which  is  the  Fourier  spectmm 
E'(A:^,z),  versus  /(to/c).  In  Figure  5,  most  Fourier  components  are  concentrated  about 
=  0.  However,  for  the  thick  curve  (z  =  0)  there  are  two  peaks  about  k^  l{(olc)  =  ±1.1 
and  we  interpret  these  as  wave  numbers  associated  with  the  surface  plasmon.  For  the  dot 
curve  (z  =  -500A.),  we  see  there  are  no  similar  peaks. 


Normalized  Wave  Number,  k  /(co/c) 


FIGURE  5.  Electric  Field  Intensity  Versus  Wave  Number 
k^  for  oo/c  =  1.5  Inverse  Microns  and  p-Polarization.  The 

solid  and  dot  curves  are  for  z  =  0  and  -500X,  respectively. 
The  solid  curve  shows  two  sharp  peaks  centered  about 

1<±A:^  <1.2  which  indicate  coupling  with  surface  plasmons 
which  propagate  along  the  jc  direction.  The  remaining 
Fourier  spectrum  is  associated  with  propagation  of  guided 
waves  in  the  planar  waveguide  along  the  z  direction. 
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x-Coordinate 


FIGURE  6.  Electric  Field  Intensity  Versus  x  Coordinate  for  s- 
Polarization  and  oVc  =1.5  Inverse  Microns.  For  reference,  the 
z  =  0  Gaussian  beam  envelope  is  indicated.  The  remaining  two 
intensity  curves  are  for  the  transmitted  field  (Equation  12)  at 

z  =  0  (thick  line)  and  z  =  -500X  (thin  line).  The  planar 
waveguide  boundaries  are  at  =  ±  1.5  micrometers  as  shown 
by  the  vertical  dash  lines. 


The  result  for  s-polarization  is  shown  in  Figure  6.  Again,  the  Gaussian  beam  at  z  =  0 
is  shown  for  reference.  The  vertical  dashed  lines  indicate  the  boundaries  of  the  waveguide 
channel  where  |x|<1.5  micrometers.  The  thin  and  thick  curves  are  the  electric  field 

intensity  at  z  =  0  and  z  =  -500A,,  respectively.  Again,  it  is  seen  that  there  is  significant 
coupling  into  the  guided  wave  channel  where  the  reflectance  and  transmittance  are  0.282 
and  0.718,  respectively.  Note  that  for  this  polarization,  there  are  no  evanescent  surface 
modes  and  this  solution  propagates  in  both  the  x  and  z  directions.  The  electric  field 
intensity  is  peaked  about  x  =  0. 


CONCLUSIONS 


We  have  described  a  method  of  calculating  the  reflected  and  transmitted  fields 
associated  with  a  Gaussian  intensity  beam  incident  endwise  upon  a  planar  waveguide 
structure.  The  numerical  results  indicate  that  a  majority  of  incident  energy  can  be  coupled 
into  the  planar  structure.  The  results  account  for  guided  modes  which  propagate  down  the 
planar  guide  and  for  surface  modes  associated  with  the  metal- vacuum  interface. 
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